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ABSTRACT 

Context. With the arrival of the next generation of ground-based imaging interferometers combining from 4 to possibly 6 telescopes 
simultaneously, there is also a strong need for a new generation of fringe trackers able to cophase such arrays. These instruments have 
to be very sensitive and to provide robust operations in quickly varying observational conditions. 

Aims. We aim at defining the optimal characteristics of fringe sensor concepts operating with 4 or 6 telescopes. The current detector 
limitations impose us to consider solutions based on co-axial pairwise combination schemes. 

Methods. We independently study several aspects of the fringe sensing process: 1) how to measure the phase and the group delay, 
and 2) how to combine the telescopes in order to ensure a precise and robust fringe tracking in real conditions. Thanks to analytical 
developments and numerical simulations, we define the optimal fringe-sensor concepts and compute the expected performance of the 
4-telescope one with our dedicated end-to-end simulation tool sim2GFT. 

Results. We first show that measuring the phase and the group delay by obtaining the data in several steps (i.e. by temporally 
modulating the optical path difference) is extremely sensitive to atmospheric turbulence and therefore conclude that it is better to 
obtain the fringe position with a set of data obtained simultaneously. Subsequently, we show that among all co-axial pairwise schemes, 
moderately redundant concepts increase the sensitivity as well as the robustness in various atmospheric or observing conditions. 
Merging all these results, end-to-end simulations show that our 4-telescope fringe sensor concept is able to track fringes at least 90% 
of the time up to limiting magnitudes of 7.5 and 9.5 for the 1.8- and 8.2-meter VLTI telescopes respectively. 



Key words. Techniques: high angular resolution - Techniques: interferometric 
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1. Introduction 

The sensitivity of ground-based interferometers is highly lim- 
ited by the atmospheric turbulence and in particular by the ran- 
dom optical path difference (OPD) between the telescopes, the 
so-called piston. By making the fringes randomly move on the 
detector, the piston blurs the interferometric signal and pre- 
vents from using integration times longer than the coherence 
time of the atmosphere tq (typically a few 10 ms in the near in- 
frared). To reach their ultimate performance and increase their 
number of potential targets, interferometers need fringe track- 
ers, i.e. instruments dedicated to measuring and compensating 
in real-time the random piston. By keeping the fringes locked 
with a precision better than Xj 10, they ensure a fringe visibil- 
ity loss lower than 20% with integration times of a few sec- 
onds. Up to now, fringe trackers had to cophase array up to 3 
telescopes by co mbining 2 base l ines (e.g., the FINITO fringe 
tracker at VLTI; iGai et al.ll200l iLe Bouquin et al.ll2009l) . The 
new genera tion of interferometri c instru ments, such as MIRC 
at CHARA dMonnier et al 1l2004[). MR OI dJurgenson et alj2008l) 
or GR AVITY dGillessen et al l 1201 Oh . MATISSE dLopez et al] 
120081) and VSI dMalbet et all 120081) at the VLTI, requires to 
cophase arrays of 4 and possibly 6 telescopes, raising new fringe 
tracking challenges. This paper aims at defining the optimal con- 
cept of fringe sensor for such arrays. 

This study is focused on solutions based on co-axial pairwise 
combination of the light beams, as currently used in existing and 
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planed fringe-tracker such as FINITO, CHAMP and GRAVITY. 
The reason is that fringe sensing is generally carried out in the 
detector-noise limited regime and that multi-axial combination 
requires a larger number of pixels than pairwise co-axial com- 
bination. Additionally, we consider only the concepts providing 
measurements of both the phase delay (phase of the interfero- 
metric fringes) and the group delay (position of the white-light 
fringe). Indeed, the group delay resolves the 27r ambiguity on the 
phase and is mandatory to ensure an efficient and robust fringe 
tracking. 



To define the optimal 4- and 6-telescope fringe sensor con- 
cepts based on the co-axial pairwise combination, we study 3 
independent points. In Section [2] we study the phase estima- 
tor. We compare two different implementations of the ABCD 
fringe coding depending on whether the ABCD samples are ob- 
tained simultaneously or sequentially. In Section [3] we study the 
two possible ways to measure the group delay, either by tempo- 
rally modulating the OPD or by spectrally dispersing the fringes. 
In Section [4] we compare the efficiency of beam combination 
schemes with various degrees of redundancy (that is forming all 
the possible baselines of the array or not). We show that the re- 
sult is a tradeoff between precision and operational robustness. 
Finally in Section [5] we merge the results of the 3 previous sec- 
tions to define the optimal concept in the 4-telescope case. We 
perform a detailled estimate of its performance in the VLTI en- 
vironnement. 
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Fig. 1. The ABCD estimator. Left: conceptual representation of 
the 4 phase states sampling the fringes. Right: the measured 
phase states functions of the time for static (top) and temporally 
modulated (bottom) ABCD. The total integration time is to. 

2. Phase estimation 

Measuring the phase is essential for a fringe tracker in order to 
stabilize the fringes and to cophase the array within a fraction 
of wavelength. In this section, we therefore consider we are in 
a cophasing/phase tracking regime in which the group delay is 
known. We compare the precision of two different implemen- 
tations of a phase estimator depending on whether the required 
measurements are simultaneous or not. The simplest and most 
efficient w ay to measure the fringe phase is the so-called ABCD 
estimator (Sh ao et ai1ll988l) . It consists in sampling 4 points in 
quadrature in the same fringe (see Fig. [TJ left), so that the real 
and imaginary parts of the coherent signal are extracted: 



A - C cc V cos (p 
D-B a Vsin <p 



(1) 



where V and are the fringe visibility and phase respectively, 
the cotangent of the latter being then estimated by: 

D-B 



tan0 



A C 



(2) 



Considering a total integration time fo to obtain a phase esti- 
mation, there are two possible ways to perform the ABCD mea- 
surements (Fig.[TJ right): 

- Temporal ABCD: it consists in temporally m odulating 
the O PD like in the cases of FINITO at VLT I dGai et all 
12004 . CHAMP at CHARA dBereer et al.ll2006h or the Keck 
Interferometer fringe tracker (IColavita et al. 20101) . We will 
consider in the following an implementation using a sam- 
pling of both outputs of a beam-splitter (in phase opposi- 
tion) simultaneously. This allows the recording of two phase 
states A and C (in phase opposition) from t = to fo/2, and 
the B and D phase states by adding a temporal n/2 phase and 
recording between t = to 12 and fo. This way one can gener- 
ate an ABCD fringe coding (see Fig. [1] right and bottom). 
There is consequently a fo/2 time delay between the (A,C) 
and (B,D) samples. Other possible implementations (for in- 
stance at the Keck Interferometer fringe tracker) consider a 
continuous modulation over 1 fringe and only use one of the 
two interferometric outputs to measure the phase. Providing 
an exhaustive comparison between possible temporal algo- 
rithms is out of the scope of the paper but might lead to select 
a different implementation. 

- Static ABCD: with this method we simultaneously measure 
the four phase states from t = to fo- Thi s method is im- 
pleme nted in the PRIMA FSU at the VLTI dSahlmann et alJ 
l2009h and is expected to be used on future instruments such 
as GRAVITY. In this case, there is no time delay between the 
ABCD samples. 



In both cases the same signal-to-noise ratio (SNR) is achieved 
since the same number of photons is collected. The static ABCD 
requires to make twice more measurements simultaneously, so 
that the output flux is divided by 2, but each pixel integrates the 
signal twice as long. However the temporal and static ABCDs 
are not fully equivalent in real conditions because of atmospheric 
and/or instrumental disturbances. We now compare them by tak- 
ing into account such effects. 

2. 1 . Phase measurement errors 

When considering piston or photometric disturbances, the phase 
quadratic error cr 2 decomposes into the sum of two terms: 



cr„;„ + cr. 



del 



(3) 



The first one is the noise due to the interferometric si gnal detec- 
tion c r s j g which includes detector and photon noises (IShao et aTl 
119881) . The second one, the so-called delay noise cr de i, is due 
to external disturbances (piston or photometric variations) that 
combine with a delay between the ABCD measurements. By 
definition, the temporal ABCD is affected by such a noise, but 
not the static ABCD, since the four measurements are simul- 
taneous. As this noise is an additionnal term, independent of 
the source brightness, we can already anticipate that it limits the 
phase measurements precision at high flux. 

2.1 .1 . Detection noise 

While integrating the signal, the fringes move slightly because 
of the atmospheric piston. Their contrast is attenuated by a fac- 
tor exp (cr 2 ((f) p , fi)/2), where <p p is the piston phase and cr 2 (<p p , ti) 
its variance for an integration time t\ . The integration time per 
phase state is twice larger in the static case than in the tempo- 
ral case (see Fig. Q] right) implying a more important contrast 
loss. Combining this effect with the expressi on of the d e tectio n 
noise for an ABCD estimator derived from IShao et al 1 (119881) . 
we obtain: 

exp (0.5 cr 2 (4> p ,to/2)) 
in the temporal case 

(4) 

exp(o.5cr 2 (0 p ,f o )) 
in the static case 



< = 2 



4cr 2 e +K 
~V 2 K 2 ~ 



where cr e is the read-out noise in electrons per pixel, V is the 
fringe contrast and K is the number of photo-events for a total 
integration time fo. The left term corresponds to the sum of the 
detector and photon noises respectively. 

2.1 .2. Delay noise 

Delay noise is the consequence of piston and photometric vari- 
ations between the (A,C) and (B,D) measurements, and there- 
fore only affects the temporal method. These disturbances can 
be induced by the atmosphere (piston and scintillation) or by the 
instruments (vibratio ns). Since at mospheric piston and scintilla- 
tion are independent dFriedll 19661) . we can decompose the delay 
noise in two terms due to the piston and the scintillation respec- 
tively: 

^del = ^pist + °ici ( 5 ) 

To compute these noises we assume the disturbances are con- 
stant while integrating the (A,C) signal, and suddenly change 
while integrating (B,D). 
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Condition 


Excellent 


Good 


Medium 


Bad 


Seeing [arcsec] 


0.46 


0.55 


0.64 


1.10 


t [ms] 


8.7 


3.1 


2.7 


2.0 



Table 1. Typical seeing and atmospheric coherence time To for 
the different observing conditions considered. 



ATs 


Piston noise <x p j st 


Scintillation noise <r sc i 


t [ms] 


2 


4 8 


2 4 8 


Good 


A/92 


/1/60 A/35 


A/499 A/369 A/290 


Bad 


A/29 


A/19 A/12 


,1/101 /1/67 ^/37 


UTs 









t [ms] 
Good 
Bad 



12 4 12 4 

A/33 A/21 A/12 A/162 A/122 A/59 
A/21 A/14 A/9 A/101 A/52 A/21 

Table 2. Piston and scintillation noises computed from Eq. [7] 
and [10] They are expressed as a function of the wavelength (in 
the H band), for three different integration times. Atmospheric 
conditions are Good (G) and Bad (B). For more details, see Tab. 
IA.ll and rA~2l in AppendixlAl 



Piston noise - Because of the piston variation between the 
(A,C) and (B,D) measurements, the phase difference between 
them is not n/2 as it should. Taking the point in the middle of 
the interval fo as the reference, the measured signal is therefore: 



A-C «V cos(0 + (p p {t - t /4)) 
D-B cxV sin(0 + (p p (t + f /4)) 



(6) 



The comparison to the ideal signal in Eq.Q] shows that the esti- 
mated phase </> is biased if (p p has varied between measurements. 
When we take into account the piston statistics, this bias results 
in the following piston noise of variance: 



cr^ st =0.125<r 2 (c5^,fo/2) 



where o- 2 {6<p p ,tol2) is the variance of the difference of piston 
separated by to/2. Details of the computation can be found in 
AppendixlAl 

Scintillation noise - The fringe visibility depends on the flux 
imbalance between the two beams I\ and li of the interferometer. 
These unequal fluxes reduce the fringe visibility by a factor: 



h+h 



(8) 



Because of scintillation, I\ and I2, and therefore V sc j, change be- 
tween the (A,C) and (B,D) measurements. Still considering the 
middle of the interval fo as the reference, the measured signal is: 



A-C « V sci (t - t /4) cos 
D-B a y sci (f + f /4) sin < 



(9) 



By comparing this equation to the ideal signal (Eq. [TJ, we see 
that a single phase estimation is biased if V sc j varies, that is if I\ 
and/or I2 vary. Assuming the beams I\ and h to be independent 
and of same statistics, the scintillation noise is: 
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Fig. 2. Relative errors cr^/A of temporal (dash) and static (solid) 
ABCD phase estimators in H band as a function of the number 
of detected photo-events K. Black, blue, magenta and red curves 
represent Excellent, Good, Medium and Bad conditions respec- 
tively as defined in Tab. [TJ The plots are done in the case of the 
ATs for an integration time of 2 ms for the specific ABCD im- 
plementation considered here. Note that for the static ABCD, 
the black, blue and magenta curves are superimposed because of 
close performances. 



where x = (7i(f +fo/4)-/i(f-f /4))/7i(f) is the relative flux vari- 
ation between the (A,C) and (B,D) exposures, (x) its mean and 
<r 2 (x, to/2) its variance during a time fo/2. Note that to compute 
this noise, we consider the extreme case of a mean unbalance 
between the interferometric inputs equal to 10. Details of the 
calculations can be found in AppendixlAl 



(') 2.2. Performance comparison 



cr 2 .~ 0.04 cr 2 (x, t /2) 



(10) 



In order to put quantitative numbers on the previous results, 
we used data provided by ESO and collected at the Paranal 
Observatory in 2008. The FITS files contain the photometric flux 
and the fringe phase as measured by the FINITO fringe-tracker 
in the H-band. Data were collected at a frequency of 1 kHz for 
ATs and 2 kHz for UTs, and for various atmospheric conditions 
(see Tab.[JJ. We have computed the variance of the difference of 
piston and photometries separated by fo/2, for different values of 
fo. We have finally injected the results in Eq.[7]and[T0]to evaluate 
piston and scintillation noises, in atmospheric conditions rang- 
ing from Excellent to Bad (see Tab.[2]i. We note that whatever the 
conditions and the integration time, cr p i St is always at least twice 
larger than cr sci : when measuring the phase, the piston is there- 
fore far more harmful than the relative variations of flux - this 
is all the more true than we consider an extremely unfavorable 
case for scintillation noise, as explained in the previous section. 

We now compute the phase error cr^ in realistic conditions 
for the temporal and static ABCD methods. Fig.[2]represents the 
phase error relative to the wavelength (i.e. cr^/A) in H-band with 
both methods. It clearly shows that the static ABCD outperforms 
the temporal one in almost all regimes. It is only in the photon 
poor regime and in bad conditions that modulating the fringes is 
a little more efficient, that is when the fringe contrast attenuation 
on the static ABCD becomes important. Yet regarding the large 
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phase error (cr^ > /1/10, see Tab. 0, phase tracking would be 
very poor - if possible - in such conditions. 

In the photon rich regime, the plateau for the temporal 
method is due to the delay noise. For the 1.8-m Auxiliary 
Telescopes (ATs) at the VLTI, it has an almost null impact on 
phase tracking in good conditions (cr^ < A/35) even for in- 
tegration times as long as 8 ms. In bad conditions with inte- 
gration times longer than 2 ms there could be some limitations 
(ct-0 > A/20) depending on the actual implementation of the tem- 
poral AB CD. 

Observations on the 8-m Unit Telescopes (UTs) show 
a higher piston noi se, partly due to instrumental vibrations 
(iDi Lieto et al. 2008): in good conditions, the noise level is simi- 
lar to the one of ATs in bad conditions. Passing from good to bad 
conditions, the integration time has to be divided by 2 to main- 
tain the performance in a photon rich regime. In particular, in 
bad conditions and to > 4 ms, the noise level is higher than Aj 10 
whatever the source, and phase tracking could be hardly possible 
with a temporally modulated ABCD. This probably explains the 
difficulty of the FINITO fringe-tracker to close the loop on the 
UTs for faint objects. 

In conclusion, with a temporal phase estimator, the fringe 
tracking capabilities are compromised in bad atmospheric con- 
ditions and on faint sources requiring long integration times. 
Therefore, from a performance point of view, a static method 
should be preferred thanks to its lower sensitivity to distur- 
bances. 

3. Group delay estimation methods 

The group delay (GD) is a measurement complementary to the 
phase and is mandatory to ensure an efficient fringe tracking. 
Indeed, a phase estimator only determines the fringe position 
modulo 27r. The GD lifts this ambiguity (see Fig. [3]). It allows 
to find and recover the position of maximum contrast, therefore 
providing the highest SNR. This is of particular interest when the 
fringe-tracking is unstable and/or when unseen fringe jumps oc- 
cur regularly. Moreover, monitoring both the GD and the phase 
allows to determine t he amount of dispersion ind uced by atmo- 
spheric water vapor dMeisner &LePoolell200l . This is done 
routinely at the Keck Interferometer for cophasing i n N-band 
while measuring the phase and group delay in K-band dColavital 
I2010h . 

1(A) and V(A) being the flux and the complex visibility of the 
interferometric signal, the coherence envelope is linked to the 
complex coherent flux I(A)V(A) through a Fourier transform: 



E(x) a 







I(A)V(A), 



dA 



(11) 



where x is the OPD. Consequently it is possible to estimate the 
group delay with two different methods: 

- The temporal method estimates the GD by measuring the 
envelope amplitude (in other words the fringe contrast) E(x) 
at several points around its maximum by modulating the 
OPD. Since the phase needs to be measured at the same time 
to ensure fringe tracking, the OPD is modulated near the en- 
velope center to keep a high SNR. This method is currently 
used in FINITO and CHAMP. 

- The spectral method uses the Fourier relation between 
the coherent spectrum I(A)V{A) and the coherence envelope 
E(x). The coherence envelope is recovered by measuring the 
coherent spectrum over few spectral channels. This method 
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Fig. 3. Top: example of polychromatic fringes (solid line) with 
longitudinal dispersion, modulated by the coherence envelope 
(dashed lines). Bottom: corresponding phase and group delay 
measurements (in blue and red respectively) presented in mi- 
crons. 



has b een successfully implemented at PTI ( Colavita et al.l 
Il999h . and more recently in PRIMA (ISahlmann et al.ll2009l) 
and in the KI fringe tracker (IColavita et al.ll2010T) . 

We could not obtain a realistic analytical description of these 
group delay estimators. Therefore we decided to compare them 
with Monte-Carlo simulations taking into account atmospheric 
disturbances. 



3.1. Description of the simulations 

We want to fairly compare both methods, so that: 

- We use the same fringe coding, i.e. a static ABCD because 
of its lower sensitivity to disturbances (see the previous sec- 
tion). 

- The signal is integrated during the same amount of time so 
that each method collects the same amount of photons and is 
prone to the same disturbances. 

- In both cases, the group delay is estimated in the same way 
by fitting an envelope model to the processed data. This al- 
lows a comparison of the intrinsic quality of the data for both 
methods. There are obviously many other ways to estimate 
xqu from a set of data , but we assume that this is a second 
order problem. Indeed. IPedretti et al.l (l2004l) compared three 
different algorithms to estimate the group delay with a tem- 
poral method and noted only little differences on the perfor- 
mance, eve n with an algorithm as sophisticated as the one 
proposed bv lWilson et al l d2004l) . 



N. Blind et al.: Optimized fringe sensors for the VLTI next generation instruments 
T emporal metho d Dispersed method 



5 



Coherent 




Fig. 4. Conceptual representation of the signal processing for group delay estimation. Temporal method (left): an envelope model is 
fitted on the 3 envelope amplitude measurements to determine the group delay. Dispersed method (right): from the spectral sampling 
of the complex coherent signal, an approximated envelope is computed with a Fourier transform operation. The envelope position 
is determined by fitting an envelope model. 



These choices made, temporal and dispersed methods can also 
be optimized in order to improve their performances. Here be- 
low, we describe the characteristics of each method. 



number of photons: therefore we make three GD estimations 
with the dispersed estimator, introducing disturbances between 
each estimation, and finally average them. 



3.1.1. Temporally modulated interferogram 

Simulations have shown that the temporal estimator is strongly 
affected by atmospheric and instrumental disturbances. Their 
effect is all the more minimized than the envelope is quickly 
scanned. Our study shows that the optimal way to proceed is to 
successively measure the fringe contrast in three different points 
over a 5-fringe range (OPD equal to -2.5A, and 2.5A near 
the envelope ma ximum). This resul t is in agreement with the 
CHAMP choice jBergeret ai]|2006l) . Once the three contrasts 
are measured, they are fitted with an envelope model to deter- 
mine the group delay. A schematic overview of this method is 
displayed in Fig. [4] left. 

The input fluxes have to be monitored to compensate in real 
time for the photometric/contrast variations that occur between 
the 3 measurements. For sake of simplicity, we consider these 
photometric estimations to be noise-free. The simulated perfor- 
mance for the temporal method will thus be optimistic. 



3.1.2. Spectrally dispersed interferogram 

Thanks to the ABCD measurements, we can compute the chro- 
matic complex visibility I(A)V{A) on each spectral channel. An 
approximated coherence envelo pe is then computed by taking 
their discrete Fourier transform dColavita et al.lll999l) . It is pos- 
sible to disperse the fringes over three channels to optimize the 
sensitivity, but we decide to use five channels to enhance the 
spectral sampling and thus the robustness of the estimator (see 
Section f3T2T > - For each exposure, a set of dispersed ABCD data 
is obtained, which enables a new GD estimation. 

For a fair comparison between the temporal and the spectral 
method, they are both fed with the same disturbances and 



3.2. Linearity and dynamic range 

A reliable estimation of the group delay is of prime importance 
since it ensures the measurements to be made in the highest SNR 
area. We study in this part two quantities, the linearity and the 
dynamic range, by looking at the response xgd of both methods 
to a given OPD ramp xgd- We define the linearity r\ as the local 
slope of xgd versus xqd: 



n = 



dx GD 



(12) 



GD 



A perfectly linear estimator is such that r\ — 1 . Otherwise the 
estimator is biased and the envelope is not perfectly stabilized. 

The group delay is extremely important for the fringe track- 
ing robustness, that is the ability of the estimator to keep the 
fringes locked in the highest SNR area, in particular after a 
strong piston stroke (> 15 /urn). In practice, there are limits 
outside which the group delay estimation is highly biased and 
makes the fringe tracker diverge from its operating point. The in- 
terval between these limits corresponds to the so-called dynamic 
range (DR), which is used here to characterize the robustness of 
the estimators. In practice, the limits of the DR are reached when 
the slope of xqd versus xqd changes sign (in other words when rj 
becomes negative) or when we observe a strong wrapping effect. 

In the following paragraphs, we simulate noise-free ideal in- 
terferograms in the H-band with a sinc-shaped coherence enve- 
lope. We fit the results with two different envelope models (a 
parabola and a sine function) in order to study its impact on the 
GD estimation. The results of this study are presented in Fig. [5] 
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Fig. 5. Response of the dispersed and temporal group delay estimators (solid and dashed lines respectively) to an OPD ramp xgd 
in H band. The ideal response is represented by the large grey line. In all cases, the coherence envelope has a sine shape. Top: the 
envelope model is a sine function. Bottom: the envelope model is a parabola. Figures on the right are zoom on the central part of 
the left-hand side figures. The DR limits are represented with blue arrows on the top-left plot in the ideal case for both estimators. 



3.2.1. Temporally modulated interferogram 



the main lobdZ] and seems weakly affected by the model quality. 



In the temporal method, the envelope model is critical to ensure a 
good linearity. Using the most appropriate sine model with ideal 
interferograms (Fig. [5] top, dashed line), the linearity is excellent 
(77 = 1), but the DR is limited to 10 fringes (+8 //m), i.e. to the 
width of the central lobe. Outside this range the GD estimation is 
totally non-linear but never cross the y-axis: the fringe tracking 
loop should not diverge but it should recover the envelope center 
with difficulty, or even could risk to lock the fringes far away 
from the envelope center. 

Using a wrong envelope model (e.g., a parabola; see 
Fig. [5J bottom, dashed line) leads to a relative bias higher than 
10% (77 ~ 0.9) whatever the OPD within the dynamic range. 
Increasing the number of samples or the scan length does 
not improve the results, emphasizing that the problem comes 
from the wrong envelope model. Because of the number of 
chromatic variables (particularly the longitudinal dispersion) 
which continually vary during a night and slightly modify the 
envelope shape, the envelope model cannot be perfect and the 
temporal estimator will therefore be consistently non linear by a 
few percents. Interestingly the DR is still equal to the width of 



3.2.2. Spectrally dispersed interferogram 

On the contrary, the dispersed method is not affected by the 
envelope model (see Fig. [5] solid lines): since we sample the 
complex coherent spectrum, we can directly compute a real- 
istic coherence envelope and the fitting model has therefore a 
weak influence. Dispersing fringes on 5 spectral channels in H- 
band, the linearity is excellent (77 ~ 1) over an OPD range of 
±20yum. Beyond these points a sharp wrapping effect is observed 
(Fig- El left), marking the DR limits: the discrete sampling of the 
spectrum induces aliasing effects on the computed envelope (ob- 
tained from a discrete Fourier transform of the complex coherent 
signal, see Eq.lTTb. so that outside the DR the GD is estimated on 
a replica of the true envelope. In practice, if the GD is measured 
after such a wrap, the fringe tracker will correct the OPD in the 
wrong direction and finally lock the fringes on a point even more 



1 Simulations show that the DR can be increased with a higher num- 
ber of contrast samples and a higher scan length. But in real operations 
it also increases the influence of atmospheric disturbances, which is not 
suitable for precision purposes (see Section [OV 
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# photo-events 

Fig. 6. Relative errors ixgdM of temporal (dash) and dis- 
persed (solid) GD estimators in different atmospheric conditions. 
Black, blue, magenta and red curves represent Excellent, Good, 
Medium and Bad conditions respectively as defined in Tab. Q] 
The plots are done in the case of the ATs, for a total integration 
time of 3 ms, constituted of 3 single exposures of 1 ms. 

distant from the envelope center than previously. However, since 
we have chosen to use 5 spectral channels, the DR (±20yt/ra) is 
larger than the strongest piston fluctuations typically observed 
on a few milliseconds (~ 15/j.m). Note that working in K-band 
increases the dynamic range up to +40fim, almost cancelling 
such issues. It is actually possible to infer an expression for the 
DR with dispersed fringes. Let us assume a spectral band with an 
effective wavelength Ao and a width AA, and that the fringes are 
dispersed over channels. The dynamic range is then ideally 
(see Appendix|B]i: 

A 2 

DR = N A -± (13) 
AA 

The larger the number of spectral channels, the lower the alias- 
ing and therefore the larger the DR. This relation is in excellent 
agreement with the simulation results. 

When longitudinal dispersion is taken into account, the lin- 
earity and DR are slightly reduced because the undersampling 
of the coherent spectrum leads to a less precise envelope com- 
putation. Refining the spectral sampling with more channels im- 
proves both linearity and DR as shown by the simulations. 

In conclusion, spectrally dispersing the fringes appears to be 
the most robust method to measure the group delay. It provides 
an estimator with: 

- a good linearity without the need of a good envelope model, 
as it inherently computes a realistic envelope; 

- a large DR allowing robust operations and quick recovery 
of the fringes over an OPD range larger than typical piston 
variations. 

3.3. Group delay measurements precision 

We now compare the precision of the GD estimators as a func- 
tion of the incoming flux and of the disturbances strength. The 
simulations consist in computing noisy interferograms in H- 
band, introducing detector and photon noises as well as pis- 



ton and photometric disturbances, which are taken from actual 
FINITO data. For each simulation, we estimate a noisy GD 
(xgd)- Its statistics over several thousands of iterations gives the 
statistical error <xgd for both estimators. 

The results for ATs and an integration time of 1 ms are pre- 
sented in Fig. [6] It shows the relative error o-qd/A on the group 
delay measurements as a function of the number of photo-events 
for various atmospheric conditions. The limitation of the tem- 
poral estimator is obvious, with a plateau due to atmospheric 
disturbances (piston mainly) which acts like an independent, ad- 
ditional noise at high flux, increasing when atmospheric condi- 
tions get worse. On the contrary the dispersed estimator appears 
weakly sensitive to these disturbances. Although we have used 
favorable hypothesis for the temporal method (the required pho- 
tometric monitoring is considered noise-free), there is no regime 
in which this concept is better than the dispersed one. For UTs, 
results are similar but with stronger limitations: it appears that 
the statistical error of the temporal estimator never goes below 
A/4 with integration time as low as 1 ms whatever the conditions. 

Additionally, all the simulations show the same dependency 
of the statistical error of both GD estimators with respect to the 
incoming flux K and the visibility v|3 

9 1 

°gd a g y2 in photon noise regime (14) 

7 1 

°gd a yi in detector noise regime (15) 

Interestingly, we find the same kind of dependency than for the 
phase (Eq.|4]i in the equivalent regimes. 

In conclusion, temporally modulating the OPD to estimate 
the group delay is not competitive with the spectrally dispersed 
fringe method, both in terms of robustness and precision. This is 
in line with the conclusion of Section [2] which showed the sen- 
sitivity of temporal fringe coding to external disturbances. We 
therefore strongly conclude that a static fringe coding scheme 
dispersed across a few spectral channels should be used to mea- 
sure the fringe phase and group delay. 



4. Optimal co-axial pairwise combination schemes 

Theoretically, it is possible to cophase an array of N telescopes 
by measuring only N— \ baselines. However because of the noisy 
measurements and of the varying observing conditions during a 
night, some baselines can deliver information of poor quality, 
so that it is beneficial to have some redundancy with additional 
baselines. It is then possible to retrieve the phase on a baseline 
in several different ways, ensuring a better fringe tracking sta- 
bility. The drawback is that when the number of measured base- 
lines increases, each one is less sensitive because the flux of the 
telescopes is divided between more baselines. The sensitivity of 
the fringe sensor then depends on a competition between the in- 
formation redundancy and the sensitivity of the individual base- 
lines. The aim of this section is to determine the most efficient 
schemes with respect to their intrinsic performance and opera- 
tionnal advantages. 

Several on-going projects will work with 4 (GRAVITY, 
MATISSE) to 6 (VSI, MIRC) telescopes. Therefore we focus 
on these 2 cases, assuming that all telescopes are identical. We 
consider the following schemes illustrated in Fig. [7] 



2 These empirical relations are only valid when there is no distur- 
bance for the temporal method. 
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Fig. 7. The various conceptual schemes studied for the 4T and 6T cases (top and bottom respectively) with the associated nomen- 
clature. As explained in section |4~T1 we do not study schemes with intrinsically imbalanced photometric inputs other than the open 
ones, because of lower performance. 



- The open schemes are made up of the minimal number of 
baselines, that is iV - 1 and are noted ATX). In this case the 
interferometric outputs are intrisically imbalanced in flux, in 
order to have baselines with equivalent performance. For in- 
stance, in the 4TO case, we do not distribute 50% of the in- 
tensity of the telescope 2 onto baselines {12} and {23}, but 
~ 40% and ~ 60% respectively (see Appendix lC.2.2l for the 
details of this optimization). 

- In the redundant schemes, the flux of each telescope is 
equally divided between the same number R of baselines. 
When R — 2 the schemes are more precisely called circular. 
The nomenclature to designate them in the following is NTR, 
possibly with an additionnal letter when there are several 
possibilities for the same value of R. 

4.1. Study of the combination schemes 

We have decided to compare the various schemes on the base of 
three considerations: their intrinsic performances, their ability 
to provide the individual beam photometries without dedicated 
outputs, and their robustness to unpredictable and rapidly vary- 
ing observing conditions. 

4.1 .1 . Performance study 

The principle of our analys is is similar to the o ne led for the 
GRAVITY fringe tracker bv lHouairi et all (I2008I) . It consists in 
computing the vector of the optimal optical path estimators x 
used to drive the delay lines, from the noisy and possibly re- 
dundant phase information <f>. These quantities are linked by the 
interaction matrix M which is known: 



(16) 



With redundant schemes, the system is overdetermined so that 
we use a x 1 minimization procedure to compute the control ma- 
trix W and then x: 



(17) 



Since <p is noisy, we have to take into account the error on the 
measurement when computing W, in order to reduce the impact 



of the noisiest baselines and prevent the solution from diverg- 
ing. The quantity of interest is finally the error cry on the cor- 
rected differential pistons calculated for each baseline [ij] with 
respect to a reference noise o"o, which corresponds to the error 
of a simple t wo-telescope interf erometer. The expression of <tq is 
derived from 
iTatulli et al] I 



Shao et al. I (fl988l) or, in a more general form, from 
2010h . It depends on the considered noise regime, 



so that the detector and photon noise regimes can be indepen- 
dently studied: 



del 
"0 

phot 
^0 



A 

K~V 

B 



(18) 
(19) 



where A and B are proportionality factors depending only on 
the fringe coding, so that this study is independent on the phase 
and the group delay estimators used. Results for the different 
schemes are therefore perfectly comparable within the same 
regime. Note that the above expressions also agree with our pre- 
vious results concerning the group delay (Eqs.[T4"land[T5V 

To analyze the behavior of the different schemes in realistic 
conditions, we consider the following three cases: 

- Ideal case: all the baselines are strictly equivalent in terms 
of flux and visibility. 

- Resolved source case: one baseline of the array is highly 
resolving the source (cases e.g. of an asymmetric source or 
of a very long baseline). To study this case, we set the fringe 
visibility to 0.1 on one particular baseline, and to 1 on the 
others. 

- Low flux case: the flux of one telescope is set to one tenth of 
the others, to simulate a quick variations of flux (e.g. scintil- 
lation) or a technical problem. 

The results for these three cases are presented in Tab. to |5] 
showing the relative error ey = cry/o"o on the corrected piston 
for the various baselines {ij}. 

In the ideal case (Tab. [3]), the redundancy slightly degrades 
the performance in the detector noise regime (because the sig- 
nal is coded on a larger number of pixels) but does not im- 
pact the performance in the photon noise regime. The differ- 
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Detector noise 


Photon noise 


Scheme 


E u 


e a 


4TO 


1.6 


1.3 


4T2 


1.7 


1.2 


4T3 


2.1 


1.2 


6TO 


1.8 


1.4 


6T2 


1.8 


1.3 


6T3A 


2.2 


1.3 


6T3B 


2.2 - 2.3 


1.3 


6T4 


2.6 


1.3 


6T5 


2.9 


1.3 



Table 3. Results of the performance study in the ideal case, 
where all the baselines are equivalent and noted [ij}. 



ences are at maximum of the order of 30% between the vari- 
ous schemes. Open and circular schemes provide similar perfor- 
mance. However, in the open schemes, the flux is not divided 
equally between the various baselines to reach an optimal SNR 
(see Appendix lC.2.2b . Although the baselines at both ends of the 
array receive roughly 40% more photons than the others, they 
are affected by a photometric imbalance, leading to a fringe con- 
trast loss of roughly 10% (i.e. an SNR loss around 20%): this 
points to the fact that the input photons are not optimally used. 
On the other hand the schemes with more baselines benefit of 
some redundancy. These facts explain why open schemes are 
slightly less sensitive in the photon noise regime than redudant 
- and balanced - ones. A simil ar conclusion concer ning open 
schemes was already reached bv lHouairi et al.l (120081) in the 4T 
case. 

In the case where a baseline resolves the target (Tab. |4}, the 
benefit of redundancy clearly appears. Indeed, whereas the mea- 
surement error on the resolving baseline strongly increases with 
open schemes, the performance degradations are well contained 
with the redundant ones. There is still a significant improvement 
between R — 2 and 3, but only limited differences between more 
redundant schemes. 

When a telescope has a reduced flux (Tab. [5j, the overall 
results do not significantly vary between the various schemes. 
Having a minimal redundancy (R = 2) appears optimal in the de- 
tector noise regime, since more baselines induce a larger overall 
read-out noise. In the photon noise regime, redundant schemes 
have very close performances and are slightly more efficient than 
the open ones. Hence circular scheme should be favored with re- 
spect to open ones and the use of more redundant schemes is not 
essential from the performance point of view. 

Taking into account the relatively close performance be- 
tween the redundant concepts and regarding their instrumental 
complexity (number of baselines to be coded, optical transmis- 
sion, etc.), schemes with R = 2 or 3 should be favored. 

4.1.2. Extracting the photometry 

The knowledge of the photometry is theoretically not mandatory 
to measure the fringe phase. However, a real-time photometric 
monitoring is very useful during operation: it provides an addi- 
tional diagnosis in case of flux-related issues and it allows the 
image quality to be optimized in all beams simultaneously (oth- 
erwise the only way to optimise the flux of each telescope is 
to optimise them sequentially). Moreover, the knowledge of the 
photometries allow the fringe visibility to be computed in real- 
time, revealing possible technical issues (or even astrophysical 
"issues" such as unknown binaries). 



Detector noise Photon noise 



Scheme 


en 


e U 


en 




4TO 


16.2 


1.6 


13.1 


1.3 


4T2 


3.4 


2.0 


2.4 


1.4 


4T3 


3.0 


2.1 -2.4 


1.7 


1.2- 1.4 


6TO 


18.1 


1.8 


13.6 


1.4 


6T2 


4.3 


2.0 


3.1 


1.4 


6T3A 


3.3 


2.3 - 2.4 


1.9 


1.3 - 1.4 


6T3B 


3.2-3.6 


2.2 - 2.5 


1.8-2.1 


1.3 - 1.4 


6T4 


3.3 


2.6-2.8 


1.7 


1.3 - 1.4 


6T5 


3.5 


2.9-3.1 


1.6 


1.3 - 1.4 



Table 4. Results of the performance study when the baseline 
{12} resolves the source. The other baselines are noted [ij] and 
are roughly equivalent. 



Detector noise Photon noise 



Scheme 




e a 




e <j 


4TO 


5.1 


1.6 


3.1 


1.3 


4T2 


4.7 


1.9 


2.5 


1.4 


4T3 


5.6 


2.4 


2.5 


1.4 


6TO 


5.7 


1.8 


3.2 


1.4 


6T2 


4.8 


1.9 


2.6 


1.4 


6T3A 


5.8 


2.2 


2.5 


1.4 


6T3B 


5.7-5.8 


2.2-2.5 


2.5 - 2.6 


1.2- 1.4 


6T4 


6.6 


2.7 - 2.9 


2.5 


1.3 - 1.4 


6T5 


7.3 


3.1 


2.5 


1.4 



Table 5. Results of the performance study in the flux drop-out 
case. The pupil 1 has a low flux and the related baselines are 
noted {1 j}. The unaffected baselines are noted {ij}. 



Some of the schemes that we study allow the instantaneous 
photometry to be extracted on each pupil without the need of 
dedicated photometric outputs. We found that, in the context of 
pairwise combinations, the photometry can be recovered from 
the fringe signal itself for every pupil that is part of a closed 
(sub-)array constituted of an odd number of pupils. Otherwise 
the system linking the fringe signals to the photometries is de- 
generated. Thus, the 4T2, 6T2 and 6T3A schemes cannot extract 
the photometry since they only contain rings of 4 and/or 6 tele- 
scopes, whereas the 4T3 and 6T3B can, since there are triangu- 
lar sub-arrays. This is summarized in Tab. |6] Note that for arrays 
with an odd number of telescopes, circular schemes (R = 2) al- 
ways allow the photometry to be directly estimated. 

4.1.3. Robustness 

When observing unknown asymmetrical sources, like well re- 
solved binary stars, unpredictable baselines can exhibit very low 
visibilities, changing with a time scale of less than one hour (see 
Fig.[8]for an example). The fringe position may then become im- 
possible to measure on some baselines, leading to a possible dis- 
continuity in the array cophasing. The case of a resolved source 
previously studied (see Tab. HJ is an example of such a situa- 
tion: when one baseline highly resolves the source, the compar- 
ison between the open schemes and the redundant ones clearly 
shows the benefit of having additional baselines. If we now as- 
sume that two baselines fully resolve the source, the schemes 
with R > 3 provide better performances than open and circular 
schemes, and so on. In general, redundancy allows bootstrap- 
ping to be performed and therefore the tracking stability to be 
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Scheme 


4TO 


4T2 


4T3 


6TO 


6T2 


6T3A 


6T3B 


6T4 


6T5 


Photometries ? 


no 


no 


yes 


no 


no 


no 


yes 


yes 


yes 



Table 6. Ability of the combination schemes to provide the inputs photometries without dedicated outputs. The schematic represen- 
tation of the schemes can be found in Fig. [7] 




V [meters] Time (hours) 



Fig. 8. Predicted fringe contrast when observing a binary star with equal fluxes and a separation of about 10 mas with the four UTs 
and a fringe sensor working in the H-band. The left panel shows the (u,v) tracks overlaid on the fringe contrast from the model. 
The right panel shows the fringe contrast versus time for 4h for each baseline. The figures have been made with the aspro public 
software from JMMC. 



increased along an observation night, so that schemes with a high 
number of baselines are favored. 

4.2. Choice of the combination schemes 

The various schemes studied here provide similar performances 
in an ideal situation. When considering more realistic condi- 
tions, the benefit of the redundancy clearly appears, by im- 
proving the tracking robustness in various observing conditions. 
Additionally, among all the schemes, some provide the input 
fluxes in real time without the need of dedicated outputs, which 
is extremely useful for the state machine. We conclude that the 
best compromises between robustness and sensitivity are the 
4T3 and 6T3B schemes. Because of their similar performance 
and their easier practical implementation, we also consider that 
the circular schemes 4T2 and 6T2 are suitable, if monitoring the 
photometric fluxes is not required. In the 4 T case, these conclu- 
sions are in agreement with the results of iHouairi et al. 1 (12001 
for the dedicated fringe tracker of GRAVITY. The results in 
the 6T ca se are also in agre ement with the choices made for 
CHAMP <Bereeretal.1 120061) at the CHAR A array with a 6T2 
configuration, even though we favor a scheme with more base- 
lines for robustness purposes. 

Despite the fact we study only 2 cases (4 and 6 pupils), it 
appears to be the trend that, in the context of pairwise combi- 
nations with an even number of telescopes, an optimal fringe 
sensor should measure either N or 3N/2 baselines {R-2 or 3 
respectively) depending on the need for photometries. With an 
odd number of telescopes, circular schemes should be optimal 
thanks to their capability to directly monitor the photometry. 

5. Estimated performance of the chosen concepts 

Now that the optimal fringe sensing concepts have been iden- 
tified, we study their on-sky performance within the VLTI in- 



frastructure. To this aim, we have developed a dedicated soft- 
ware simulation tool called Sim2GFT (2GFT standing for the 
"2nd Generation Fringe Tracker" of the VLTI). This simulator, 
consisting in a set of IDL routines, aims at performing realistic 
simulations of future observations with the 2GFT fringe sensor 
and to evaluate its performance in terms of residual piston jitter 
after closed-loop control. In the rest of this section, we assume 
that singl e-mode fibers are u sed to filter the input wavefronts, 
following lTatulli et all <l2010t) . 

5.1. The Sim2GFT simulator 

Sim2GFT is lar gely based on the GENIEsim software 
(lAbsil et al.H2006h . and therefore follows the same architecture 
and philosophy. The simulations are taking into account all ma- 
jor contributors to the final performance, from the atmosphere 
and the telescopes down to the fringe sensor and delay lines. 
The signal-to-noise ratio on the phase measurement in the fringe 
sensor is mainly driven by the amount of coherent and incoher- 
ent photons (including the atmospheric and instrumental ther- 
mal emission), and by the way they are distributed on the de- 
tector. In order to properly estimate the amount of coherent and 
incoherent photons, all the VLTI and 2GFT subsystems are de- 
scribed by their influence on the intensity, piston, and wavefront 
quality of the light beams collected by each telescope. The esti- 
mated instrumental visibility within the fringe sensor takes into 
account the visibility loss due to piston jitter, atmospheric refrac- 
tion, intensity mismatch between the beams due to atmospheric 
turbulence (scintillation), and longitudinal dispersion in the de- 
lay lines. In the case of piston jitter, a semi-empirical law based 
on on-sky FINITO data is used to include both the effect of at- 
mospheric piston and vibration-induced piston. Another key el- 
ement in the simulation is the coupling of the light beams into 
single-mode fibers, which we estimate by separating the contri- 
bution of tip-tilt (through the overlap integral between an offset 
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Airy pattern and the fiber mode) and higher order aberrations 
(through the estimated Strehl ratio — without tip-tilt — that acts 
as a multiplicative factor). 

The operation of 2GFT is closely related to the detector read- 
out scheme. Assuming a HAWAII-2RG focal plane array, we 
consider that the ABCD outputs of all baselines are spread on 
a single detector line, and that the spectral dispersion is per- 
formed on five contiguous detector lines. The detector is read 
line by line, with a read-out time that depends on the particu- 
lar arrangement of the ABCD outputs on the lines (it amounts 
to 201 fis for our design). Deriving a reliable estimation of the 
phase and group delay requires the five spectral channel to be 
usecQ. However, it must be noted that the phase and group de- 
lay estimations can be updated each time a new detector line 
is read, although it will be partly redundant with the previous 
estimation — this corresponds to the slid ing-window estimati on 
already in use at the Keck fringe tracker (IColavita et al.ll2010T) . 

The closed-loop behaviour of the fringe tracker is simulated 
by feeding back the fringe sensor phase delay measurements to 
the VLTI delay lines, using a simple PID as a controller. Group 
delay measurements are not explicitly used in our simulations, 
although in practice they will be used to make sure that fringe 
tracking is performed on the appropriate (white-light) fringe. 
The closed-loop simulation relies on a frequency-domain de- 
scription of the input disturbance (by its power spectral density) 
and of the subsystems (by their transfer function). The repeti- 
tion frequency of the loop and the controller gain are optimised 
as a function of the input photon flux and atmospheric piston to 
produce the smallest possible piston residual at the output of the 
closed loop. In order to ensure a stable fringe tracking, we re- 
quire the sensing noise to be smaller than 100 nm RMS for 90% 
of the measurements on any individual baseline, which would 
correspond to an SNR > 4 on the fringes in K band for 90% of 
the measurements. 

In the following sections, we describe the estimated per- 
formance for fringe sensing and fringe tracking of the 4T3 re- 
dundant concept with ABCD encoding on five spectral chan- 
nels over the K band (from 1.9 to 2.4 /mi). The estimations are 
based on an expected K-band transmission of 3% for the whole 
VLTI/2GFT instrument. 

5.2. Fringe sensing performance 

End-to-end simulations of VLTI/2GFT have been performed us- 
ing the 1.8-m Auxiliary Telescopes (ATs) for a K0III star lo- 
cated at various distances ranging from about 10 pc to 2kpc, 
in standard atmospheric conditions: seeing s = 0.85", coher- 
ence time to = 3 ms, outer scale £. out = 25 m, and sky tem- 
perature r s k y = 285 K. The target star is assumed to be located 
close to zenith. For each magnitude, the closed-loop repetition 
frequency has been chosen as high as possible within hardware 
limitations (< 4kHz), while keeping the average fringe sensing 
noise smaller than 100 nm RMS on all measured baselines. 

Fig.[9]illustrates the sensing noise per baseline as a function 
of stellar magnitude (black diamonds). The respective contribu- 
tions of photon noise and detector noise are represented by dot- 
ted and dashed lines. On the bright-side end of the plot, photon 
noise dominates the noise budget. The increase in photon noise 
from K = 1.5 to K — -2 is due to the star being (strongly) re- 



3 To perform a phase delay estimation with the ABCD scheme, one 
spectral channel is theoretically sufficient. However, for a better robust- 
ness to dispersion effects, we assume that the information from all five 
spectral channels is needed and will be used in practice. 
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Fig. 9. Fringe sensing noise plotted at percentile 0.9 (i.e., the 
noise is actually smaller than the plotted curves for 90% of the 
occurrences) as a function of the target's K magnitude (or of the 
stellar flux in photons per second at the entrance of 2GFT) in the 
4T3-ABCD case, assuming a K0 III target and using the A0-G1- 
K0-I1 quadruplet of ATs at the VLTI. The fringe tracking loop 
is operated at its maximum frequency as long as the fringe sens- 
ing noise per baseline remains < lOOnm RMS for 90% of the 
measurements on any individual baseline. The closed-loop repe- 
tition frequency is reduced to maintain this level of performance 
otherwise (this happens for K > 5 in the present case, as also 
shown in Fig. [Tol l, until this level cannot be reached any more 
(beyond K = 7.5 in the present case). Note that the increase in 
sensing noise for bright targets is due to the stellar photosphere 
being resolved, which reduces the available coherent flux. 



solved, which reduces the available coherent flux. Detector noise 
becomes larger than photon noise around K = 3, and the fringe 
sensing noise reaches its allowed limit (< 100 nm RMS for 90% 
of the measurements) around K — 5. For fainter magnitudes, 
Sim2GFT makes sure that the fringe sensing noise remains at 
the same level by reducing the closed-loop repetition frequency 
(i.e., increasing the integration time on the fringe sensing detec- 
tor). This is possible only until magnitude K — 7.5 in the present 
case, where a phase sensing noise of lOOnm per baseline can- 
not be reached any more for any integration time, because of the 
strong fringe blurring that appears at long DITs. The points plot- 
ted in the figure at K > 7.5 do not comply with our requirements 
any more, and have been computed for the repetition frequency 
that minimizes the fringe sensing noise (~ 33 Hz in the consid- 
ered cases). 

Also represented in Fig.|9]is the fringe sensing noise per tele- 
scope, which results from the optimized estimation of individual 
telescope pistons from all measured baselines, as explained in 
Section 14.11 The fringe sensing noise per telescope is signifi- 
cantly smaller than the measurement noise on each individual 
baseline because the estimation of the former is based on the 
information collected by multiple baselines. 

The same kind of performance study has been carried out 
in the case of the Unit Telescopes, showing a similar general 
behaviour as in the case of ATs. The only differences are: 

- the magnitude where stable closed-loop fringe tracking be- 
comes impossible, which is now around K = 9.5, 
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Fig. 10. Left: Closed-loop repetition time and time delay in the loop as a function of target K magnitude. For stars fainter than 
K = 6, the loop repetition time is increased (i.e., its frequency decreased) to ensure a sufficient SNR on the detected fringes in each 
individual measurement (until the specified SNR cannot be reached any more whatever the integration time). Right: Noise residuals 
at the output of the closed loop, for the three main contributors: fringe sensing (FSU), delay line (DL) and atmospheric noises. 



- the decrease in the coupling efficiency for stars fainter than 
V = 10, which is due to the reduced performance of the 
MACAO adaptive optics system. 

The latter effect, which is almost nonexistent in the case of ATs 
(equipped with STRAP for tip-tilt control), speeds up the drop 
of closed-loop performance at faint magnitudes. The maximum 
loop repetition frequency (~ 4 kHz) can actually be maintained 
until K =i 8.5 in the case of UTs. The presence of telescope 
vibrations in the case of UTs is taken into account in a semi- 
empirical way in our simulations, through an estimation of the 
visibility loss due to vibration-induced piston jitter, so that the 
SNR in the fringe sensing process is estimated in a realistic way. 
However, let us note that telescope vibrations are expected to 
strongly affect the residual piston jitter at the output of the closed 
loop (an effect not simulated in Sim2GFT), so that the results 
presented in right-hand side plot of Fig. [I0](in the case of ATs) 
would be significantly degraded in the case of UTs. 

5.3. Fringe tracking performance 

Fig.[l0]shows the characteristic times of the closed loop and the 
noise residuals at the output of the fringe tracking loop. The left- 
hand side plot shows that for magnitudes brighter than K — 5, 
the loop can be operated at its maximum repetition frequency 
(3.6 kHz in this case). For fainter targets, the repetition time is 
gradually increased to keep a sufficient SNR on each individ- 
ual fringe measurement. The sudden increase in repetition time 
around K — 6 is due to a modification in the loop behavior: for 
bright stars, only one spectral channel is read for each repetition 
time and the information at other wavelengths is taken from pre- 
vious repetition times, while for fainter stars all spectral channels 
are read during each repetition time (the main goal of this being 
to keep the time delajy in the loop reasonably short even at low 
repetition frequencies). The time delay is longer than the repe- 

4 The time delay of the loop is denned as the amount of time be- 
tween the middle of the overall integration time used for a phase es- 
timation (i.e., including the contribution of all spectral channels), and 
the moment when the detector read-out sequence is completed for the 
considered spectral channel. 



tition time in the bright target case, because only one spectral 
channel is read per repetition time, while the phase estimation 
uses the phase information from all five spectral channels. 

The left-hand side plot of Fig. [10] can be used to derive a 
limiting magnitude for the chosen fringe sensing concept. One 
just needs to define a repetition time threshold above which 
fringe tracking becomes inefficient. Here, we assume a maxi- 
mum allowed repetition time of 10 ms (i.e., minimum frequency 
of 100Hz)3 which gives a limiting magnitude of K = 7.5 on 
the ATs. In the case of UTs, the limiting magnitude amounts to 
K = 9.5. In both cases, this coincidentally corresponds to the 
magnitude where maintaining a phase measurement error below 
lOOnm is not possible, which indicates that a DIT of 10 ms is 
actually a sound choice to define limiting magnitudes in closed- 
loop fringe tracking operation. Note that a limiting magnitude of 
K = 7.5 in closed-loop fringe tracking with 90% locking ratio 
under standard atmospheric conditions corresponds quite well 
to what has been dem onstrated on-sky with the PRIMA fringe 
sensor unit on the ATs dSahlmann et al. 2009|). 

The right-hand side of Fig. [10] shows the noise residuals 
at the output of the fringe tracking loop, computed per tele- 
scope. Note that the fringe sensing noise residual at the out- 
put of the loop is much smaller than the actual fringe sensing 
noise (evaluated at the detection level), due to closed-loop fil- 
tering. Also note that the fringe sensing noise is always much 
smaller than the atmospheric noise under typical atmospheric 
conditions, a behavior directly related to the constraint imposed 
on the phase sensing noise per baseline in each repetition time 
(< 100 nm RMS for 90% of the measurements). For these rea- 
sons, the fringe sensing noise does not significantly affect the 
residual noise level at the output of the fringe tracking loop. The 
influence of the fringe sensor on the residual piston noise comes 
rather from its intrinsic sensitivity, which determines the maxi- 



5 For an integration time of 10 ms on the fringe sensor, the estimated 
loss of visibility due to piston jitter in standard atmospheric conditions 
is only 5% in the case of ATs, while it amounts to 28% in the case of 
UTs (an effect mainly due to vibrations). Operating at lower frequencies 
would become impractical in the case of UTs, but could be considered 
in the case of ATs (especially in good atmospheric conditions). 
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mum repetition frequency that can be reached for a given coher- 
ent flux. 

We have also performed simulations in various atmospheric 
conditions, ranging from bad (seeing eo = 1.1" and coherence 
time to = 2 ms) to excellent (eo = 0.5" and to = 10 ms). The 
influence of atmospheric conditions on the fringe tracking per- 
formance is mainly twofold: on one hand it determines the input 
atmospheric noise that needs to be corrected, and on the other 
hand it affects the amount of available coherent photons since it 
determines the injection efficiency into single-mode fibres. Our 
simulations have shown that the limiting magnitude increases 
by about 2 magnitudes between bad and excellent conditions. 
For instance, if one defines the limiting magnitude at 100 Hz, it 
varies between K = 6.2 and K = 8.5 depending on the con- 
ditions. These limiting magnitudes do not mean however that 
fringes cannot be detected at fainter magnitudes. We estimate 
that the ultimate limit for fringe detection (fringes detected for 
50% of the measurements at an SNR of 4, using a DIT of 25 ms) 
should be around K — 9.5 for ATs used in good atmospheric 
conditions. 



6. Conclusions and perspectives 

We determined the optimal 4- and 6-telescopes fringe tracker 
concepts. We showed that for realistic atmospheric conditions, 
the measurements of the various phase states (e.g., ABCD) that 
are needed to derive the fringe phase should better be done si- 
multaneously in order to limit the influence of external distur- 
bances (piston, scintillation, vibrations, etc.) on the measure- 
ment precision. Furthermore, spectrally dispersing the fringes 
allows the group delay to be evaluated with one set of contempo- 
raneous data, which (like for the phase measurement) minimizes 
the influence of disturbances. We also showed that this method 
is more robust to longitudinal dispersion effects. Therefore, we 
concluded that the optimal way to measure the fringe position 
(phase and group delay) is to perform a static ABCD fringe cod- 
ing, dispersed over about five spectral channels. 

We also demonstrated that the co-axial pairwise combination 
schemes with a moderate redundancy provide the best compro- 
mise between sensitivity and robust operations. They are less 
sensitive to varying observing conditions, and some schemes 
also allow the photometries to be directly extracted from the 
fringe signal, which is useful for the state machine. We finally 
favored the 4T3 and 6T3B schemes for 4- and 6-telescope oper- 
ations respectively. 

Merging these results, we have simulated the expected per- 
formance of the 4-telescope concept. For an efficient fringe 
tracking, with fringes locked at least 90% of the time, we ex- 
pect limiting magnitudes of 7.5 and 9.5 at K band with ATs and 
UTs respectively. These performances are close to those of sin- 
gle baseline fringe trackers currently in operation. Another im- 
portant result is that the fringe tracker ultimate performances are 
not limited by the fringe sensing measurement errors, but rather 
by the time delay between the measurement of the piston and its 
correction by the delay lines. 

Finally, in the coming years, a new generation of infrared 
detectors should be available. By providing very high acquisi- 
tion frequencies and an extremely low read-out noise at the limit 
of photon-counting, multi-axial schemes should be reconsidered 
as a possible solution for fringe-traciking, as they would not be 
limited by the large amount of pixels needed to encode the inter- 
ferometric signal. 
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Appendix A: Phase error: detection and delay 
noises expressions 

Considering an ABCD fringe coding (IColavita et al.|[l999l) . the 
phase is extracted as follows. First we have the 4 ABCD mea- 
surements in quadrature: 



A a Vcos(0) 

B a y cos(0 + tt/2) = -Vsm(q» 

C oc Vcos(4> + k) = -Vcos((/>) 

D a Vcos(4> + 3tt/2) =Vsin(0) 



(A.l) 



where V and </> are the fringe contrast and phase respectively. We 
extract the real and imaginary part of the complex fringe signal: 



\A-CoL Vcos(cf>) 
\D-B cx Vsin(0) 

and finally we estimate the phase through its cotangent: 
D-B 

tan(<j&) 



A-C 



(A.2) 



(A.3) 



We are interested here by the statistical error on the phase 
measurement, which depends on three sources of noises: detec- 
tor noise, photon noise and delay noises. Since these noises are 
statistically independent, the variance on the phase measurement 
cr? is simply the quadratic sum of these three noises: 



del 



(A.4) 



A.1. Detection noises 

The detector and photon noises terms {<r 2 det and cr 



■phot ^P^- 



tively) are derived from IShao et aT | (I1988I) for the ABCD fringe 
coding, and for sake of simplicity we put them together into the 
so-called signal detection noise a 2 -: 



del 



°~phot 



= 2 
= 2 



4^ 
V 2 K 2 
K 

Wk 2 



(A.5) 
(A.6) 
(A.7) 



where K is the number of photo-events collected during the ex- 
posure and cr 2 is the detector read-out noise. 

A.2. Delay noise 

The delay noise is due to the delay between the various measure- 
ments needed to estimate the phase and therefore only concern 
a temporal phase estimator. Because of instrumental or atmo- 
spheric disturbances (e.g. fluctuation of the differential pis ton or 
scintil lation) the phase estimation can highly biased. Since iFriedl 
(fl966h has shown that atmospheric piston and scintillation are 
uncorrected, we can study both effects independently: 



2 _ 2 2 
a del ~ Vpist + °7c, 



A.2.1 . Piston noise: cr 



(A.8) 



pisl 



We note here <p p (t) the piston term introduced by the atmosphere 
at a moment t and consider that each (A,C) and (B,D) measure- 
ment last half the total integration time to- Taking the point in 



the middle of the interval to as the reference, the interferometric 
signal writes: 



A-C ct y cos(0 + 4> p {t - to I A)) 
B - D a V sin(0 + <p p (t + k) I '4)) 



(A.9) 
(A. 10) 



We note i 



<p p (t + to/4) - <p p (t - to/4) the piston fluctuation 



between both measurements: 

A-C a y cos(0 - 64> p /2) 
B-D cx y sin(0 + 6</> p /2) 

6<p p being unknown, the phase estimator ^ is: 

B-D sin(0 + 64> p /2) 



tanc 



A-C cos(0 - 5</> p /2) 



(A.ll) 
(A. 12) 



(A.13) 



As soon as 8<p p is non null, the phase measurement is biased. If 
we consider the statistic variations of the piston, this bias can be 
considered as an additional noise. We now calculate the standard 
deviation of this phase measurement linked to the piston varia- 
tions between 2 exposures separated by a time to/2. The standard 
deviation of the piston for this time will be noted cr(5<pp, to/2). 
Assuming that the piston variations are small {o-{6<t> n , to/2) <sz 
1 rad) and using the second order expansion formula of lPapouiisI 
(Il984h . the measured phase variance writes as: 



954> p 



a- 2 (5</>p,to/2) 



(A. 14) 



<<%> 



where (5<p p ) is the mean piston variation during to/2. One shows 
then that: 



1 



cos(2 (p) 



d6c/) p 2 cos 2 (0 ■ 
Assuming that {5</> 

cr 2 0) = i cos 2 (2 , 



J2) 



sin(0 + 5<p p /2) 
cos(0 - S(p p /2) 



2\ 



(A.15) 



,) = 0, we obtain the scintillation noise: 



><T 2 (6(f) p ,to/2) 



(A. 16) 



This result depends on the mean phase position. Assuming that 
tp is uniformly distributed over [0, 27r], one finally obtains: 



'2\2 



0.125 (T 2 (54>,„tol2) 



(A. 17) 



This deviation is evaluated here by means of VLTI/FINITO data, 
and results are presented in Table IA.1I for typical integration 
times from 2 to 8 ms for ATs and from 1 to 4 ms for UTs. 



A.2.2. Scintillation noise: cr sci 

The influence of scintillation (i.e., photometric variations) be- 
tween (A,C) and (B,D) measurements is to induce fringe contrast 
fluctuations, which can bias the phase measurement. This effect 
will be studied in the same manner than in the previous section. 
Considering an ideal interferogram, the real and imaginary parts 
of the coherent signal write: 



A-C a V sci {t - t / 4) COS( 

B-D a V sci (t + t / 4) sin q 



(A. 18) 
(A. 19) 



where V SC i is the contrast attenuation term due to the photometric 
imbalance between the two beams I\ and l^. 



V sci = 



h+h 



(A.20) 
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Table A.l. Piston noise calculated with different sets of data on 
VLTI telescopes in H-band. The noise is written respectivily to 
the wavelength, for 3 different integration times. The values cor- 
respond to the worst case (<x(</>) = cr((/) p ,to/2)). Atmospheric 
conditions are: Excellent (E), Good (G), Medium (M), Bad (B). 
The corresponding observing conditions can be found in Tab. Q] 



Noting the flux variation 6, = U(t + to/ 4) - 7,(f - fo/4), the phase 
estimator writes: 



- B-D 
tan d> = = a tan < 



A-C 



(A.21) 



where: 



h+di/2 h + 62/2 h+I 2 - d 1 /2-5 2 /2 

a - -1 x (A. 22) 

\/i-<5i/2 I2-82I2 h+I 2 + d 1 /2 + 6 2 /2 

Simplifying the first and second terms by IJ2 and I\ + h respec- 
tively : 



(1 



x x /2) (1+X2/2) 1 

X 



■yi/2-y 2 /2 



<(l- Xl /2)(l-x 2 /2) l+ yi /2+y 2 /2 
with 

x t = 5i/Ij 

yt = Si/ih + h) 



(A.23) 

(A.24) 
(A.25) 



If the flux varies between both quadratures, a + 1 and the phase 
estimation is biased. If we consider the statistic variations of 
the both photometries, this bias can be considered as an addi- 
tional noise. We therefore calculate the measured photometric 
variance functions of the variance of the relative photometries 
<x 2 (x,, fo/2) between two exposures distant of to/2. We assume 
that the two pupils are sufficiently distant to be considered as 
uncorrelated, which is the case if the baseline is longer than the 
atmospheric outer scale (typically 20 m). Since the atmosphere 
follows the same statistics on both, it implies (x\) = {x 2 } = {x} 
and cr(xi, to/2) = a{x 2 , to/2) = cr(x, to/2): 



dxi 



<*i>,te> 
2 



1J7 



cr (>2,fo/2) 



cr 2 {x,tol2) 



(A.26) 



with: 



dxi 



d (a tan 0) 



1 



<*1>,<*2> 



Ox, 1 + (a tan </>) 2 



(A.27) 



In order to obtain an analytical expression of this quantity, we 
assume the flux variaitons to be small: 6j <K We do a first 
order expansion of a and only conserve the terms of the first 
order: 



a ~ (l+x 1 /4) 2 (l+X2/4) 2 (l- 3 ; 1 /2- } ;2/2) 2 
~ (l+^/2)(l + ^/2)(l- yi -y 2 ) 
~ (1 + Xl /2 + x 2 /2)(l -yi-y 2 ) 
~ 1 + Xi/2 + x 2 /2 - yx - y 2 

We expand the latter expression to simplify it: 

61 + 6 2 



1 + — + — 

+ 2/i + 2/ 2 



h +h 



And finally : 



1 + 
1 + 

1 + 



5 x l\ + b 2 l\ - {S\ + 5 2 )hh 

2(/i + / 2 )/i/ 2 
h -h 5\I 2 - 6 2 l\ 
h+h 2hh 



I /2-/1 
2 h + h 



(x\ - x 2 ) 



(A.28) 
(A.29) 
(A.30) 
(A.31) 



(A.32) 

(A.33) 
(A.34) 

(A.35) 



Therefore, considering a ~ 1 in the second term of the expres- 
sion rA23 



Noting that: 
h-h " 



I /2-/1 
2 h+h 



■ cos sin < 



(A.36) 



/2 + /1 



1-4- 



hl 



\J2 



(/l + / 2 ) 2 



1 - V 



(A.37) 



we finally show the variance of the phase measurement due to 
the photometric noise is: 

cr 2 (^) = - (sin cos 0) 2 (1 - V 2 sci ((h), (I 2 ))) cr 2 (x, t /2) (A.38) 

Note that the result depends on the mean value of the scintillating 
visibility V JCJ . Hence a perfectly balanced system should present 
a null photometric noise. This is an unrealistic effect due to our 
symetric modeling of the photometric variation with a step. In 
practice, the quick variations of photometries (i.e. during the in- 
tegration) induce a noise even for a perfectly symetric combiner. 
To obtain a more realistic value, we can consider a (worst) case 
with a mean imbalance between fluxes of a factor of 10, so that 
V sci ~ 0.57 and 1 - V 2 ci ~ 0.67. 

If we finally average this result over every realisation of <p 
(still assuming its statistics to be uniform between and 27r): 



cr 2 (4>) = 0.04cr 2 (x,f /2) 



(A.39) 



Similarly to the piston noise, we present in Table IA.2l the results 
obtained from ESO data on ATs and UTs, for different integra- 
tion times. 



16 



N. Blind et al.: Optimized fringe sensors for the VLTI next generation instruments 



ATs 



to [ms] 


2 


4 


8 




A/nyy 




i /onn 


A A 

M 


A/j4y 


i iir\ 1 

vi/JlM 


1 / 1 

vi/ lOj 


M 


Af lyo 




3/1 in 
a/ 13U 


A /I 

M 


i //inn 


1 /T77 


1/1 m 

Aj iyl 


u 

D 




A/O / 


1 /Q7 
Aj i 1 






UTs 




h [ms] 


i 


2 


4 


G 


A/162 


/I/ 122 


A/59 


M 


/I/ 107 


A/76 


A/35 


B 


A/101 


A/52 


A/21 



Table A.2. The photometric noise written respectivily to the 
wavelength in H-band, for 3 differents integration times. The 
values correspond to the worst case as defined in Eq. IA.39I 
Atmospheric conditions are: Exceptionnal (E), Good (G), 
Medium (M), Bad (B). The corresponding observing conditions 
can be found in Tab. [TJ 



where T* is the complex conjugate of T ■ Since each spectral 
channel has the same width, cr^ - cr\ — 6cr(k - I) and we finally 
get: 

E 2 (x) ^I 2 V 2 YjYj e- ilnSLT(x - XGD){k - l) (B.7) 
£=1 z=i 

The group delay is obtained when this quantity is maximum, that 
is when all the phasors in the double summation are in phase. In 
the present simple case, it is obvious it happens when x - xcd, 
which leads to: 

V(/U), e - i2KSa ' < - x - XGD) = 1 (B.8) 
And solving this equation finally gives: 

x = x GD [1 /da] (B.9) 



Appendix B: Theoretical dynamic range for the 
group delay estimation with dispersed fringes 

We analyze here the case of a dispersed estimator for the group 
delay, similar to what is implemented on PRIMA, PTI or KI. 
We remind that the coherence envelope E(x) corresponds to the 
Fourier transform modulus of the coherent signal: 



E(x) 



f 

Jo 



I(A)V(A)e i2 " XGD/A e ~ i2 * x/A dA 



(B.l) 



where x is the OPD, xqd the position of the envelope center, and 
1(A) and V(A) the source intensity and visibility, both depend- 
ing of the wavelength A. We consider a spectral band centered 
around Aq and of width AA, so that the coherence length L c of 
the wide-band interferogram is L c = A2/AA. The fringes are dis- 
persed over Na spectral channels of equal width 5 A = AA/N,\- 
In term of wavenumber, the wide- and narrow-band widths write 
Act = l/L c and 6cr = Ao-/N A . 

For sake of simplicity we consider here an ideal case, that 
is all the considered quantities are achromatic, in particular the 
source flux / and visibility V do not depend on the wavelength. 
We assume we dispose of a fringe coding (ABCD for instance) 
allowing the complex fringe signal to be computed in each 
channel k, this latter being defined as: 



Z k = I k V k e 



ilncr k x aD 



= IVe 



i2mr k x GD 



(B.2) 



where cr^ = l/A^ is the effective wavenumber on each spectral 
channel. The discrete Fourier transform of this coherent signal 
is then: 



k=\ 

Na 

— i y e -i2mr t (x-x GD ) 



k = \ 

and we finally compute the squared coherence envelope: 

E z (x) a \T(x)\ 2 

= T(x)T\x) 

k=\ I=\ 



(B.3) 



(B.4) 
(B.5) 

(B.6) 



where [ ] is the modulo symbol. In other word, by dispersing 
the fringes, we find the group delay with an ambiguity equal 
to 1 /6cr. From the definition of 5cr, it finally corresponds to an 

A 2 

ambiguity (or a dynamic range) equal to Na — ■ 

AA 



Appendix C: Noise propagation on pairwise 
combination schemes 

The study conducted in Section|4]aims at comparing various co- 
axial pairwise combination schemes looking at the phase and 
group delay measurements precision in various configurations. 
This study is based on analytical descriptions of measurement 
noises. We describe here various points which have been neces- 
sary to carry out this study but which are not essential for the 
comprehension of the results. 

C.1. Reference noise 



Thanks to lShao etall dl988h . lTatulfi et all d2010h and our study 
(Eqs.[T4land[T5Tl, we know the analytical expression of the phase 
and group delay noises, in detector and photon noise regimes and 
for co-axial pairwise combinations. They express as: 



phot 



A 

~KV 
B 



(C.l) 
(C2) 



K and V being the number of photo-events and the fringe visibil- 
ity. A and B are proportionality factors depending on the fringe 
coding, which have no influence in the following. These expres- 
sions correspond to the noise for a two-telescope (one baseline) 
instrument and are considered as noise references in the follow- 
ing. 

C.2. Individual baseline noise 

When we consider an interferometric array with more than 2 
telescopes, the flux of each telescope is distributed between sev- 
eral different baselines, increasing the noise on each baselines. 
We consider two cases here: the open and redundant schemes. 
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C.2.1. Redundant schemes 

The most simple cases are the redundant schemes in which the 
flux of each pupil is divided between the same number R of base- 
lines. Compared to a two-telescope instrument, the total flux K 
on each baseline is divided by R, so that the measurement noise 
is: 

a*" =A — 
KV 

del 



Rcr] 



-phot 



B 



(C3) 



(C4) 



We are therefore able to compare the different schemes on de- 
tector and photon noise regimes on the base of a reference noise. 




Fig. C.l. Open schemes we consider in the 4 and 6 telescopes 
cases. The nomenclature for the flux split ratio 8/ are represented 
on the figures. 



C.2.2. Open schemes 

The open schemes use the minimal number of baselines enabling 
the array to be cophased, that is N - 1 baselines. In this case the 
array is not symmetric, so that splitting the flux of intermedi- 
ate pupils into equal parts (i.e. taking 50% of their flux for each 
baseline) implies unequal performances for the different base- 
lines. In this study we want the open schemes to have intrin- 
sically equivalent baselines, that is with the same SNR on the 
fringe position measurements. To do so, we have to consider in- 
trinsically imbalanced photometric inputs for each baselines and 
we evaluate the optimal fraction of the flux to inject in the dif- 
ferent baselines. 

Considering two identical telescopes i and j, we combine 
their light by taking a fraction 8, and 8j of the incoming fluxes on 
each telescope respectively. In this case, the total flux available 
on the baseline is: 



Our goal is to have equivalent baselines, i.e., we want cry to be 
equal on the three baselines. Solving this system in detector and 
photon noise regimes leads to: 



8 de < = 0.38 



6 phot _ _ 42 



ajf = 1.62crf cr^ ' = 1.31 erf* 



6TO case For symmetry reasons we have 8\ = 8& = 1, 8% = 85 
and 83 = 84. The measurement errors on the 5 baselines write: 
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K' = K(6i + 6j)/2 



(C.5) 



and the fringe contrast V is possibly reduced because of the pho- 
tometric imbalance: 



V = V- 



5i + 5j 



(C6) 



Now considering the noise expressions in Eq. IC.ll and IC.2I we 
can easily write the measurement noises in this case, still as a 
function of our reference noises: 
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6i + 6j 
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(C7) 



(C8) 



The open schemes with 4 and 6 telescopes are presented on 
Fig. IC.ll with the associated nomenclature in term of splitting 
ratio Si. We determine in the following their values. 



cr. 



ilei 



£ 3 



phot 1 phot 



In the same way than previously, we estimate the optimal value 
of the different 5,: 



6 d 2 et = 0.31 
<5f = 0.55 



5* ot = 0.37 
6f ( " = 0.54 



ojf = 1.81 erf cr^°' 



1.36cr; 



phot 




For the 4TO and 6TO cases we note the different values of 
6i are close in detector and photon noise regimes, so that such 
schemes are practically possible. In both regime we consider the 
same values: 6 = 0.40 in the 4T case; 62 = 0.34 and £3 = 0.54 in 
the 6T case. 



C.3. Estimating the individual fringe position and final 
measurement noise 



4TO case For symmetry reasons we consider 8\ = 64 = 1 and 
62 = 83 = 8, and therefore the measurement errors on the 3 
baselines write: 



de, 
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1 phot 



We have to estimate N— \ differential pistons in order to cophase 
the interferometric array. In practice we measure B differential 
pistons (noted 0), with B > N - 1 for redundant schemes, and 
B — N - 1 for the open ones. Noting x the vector of the N - I 
optical path estimators used to drive the delay lines, the equation 
system linking cp and x is: 



= Mx 



(C9) 
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where M is the so-called interaction matrix, which is known. We 
now need to inverse this system by computing the control matrix 
W: 

x = W0 (CIO) 

For the redundant schemes, M is rectangular and we compute W 
on the base of a singular value decomposition of M. We therefore 
solve the system in the sens of a least square minimization, i.e. 
we minimize the quantity: 

X 2 = \0-Mx\ 2 (C.ll) 

However the measurements (f> are noisy and we have to weight 
them to minimize the influence of the noisiest baselines. 
Considering that the measurements have gaussian statistics and 
are statistically independent, the^ 2 writes: 



XT 



■ Mi 



(C12) 



where cr is the vector of the error on the measurement fa given 
by eq. IC.3l and lC.4l depending on the noise regime. We modify 
in consequence the differential phase vector cf> and the interaction 
matrix M as follow: 

My -» Mij/<ru j € [l,N- e [1,B] (C.13) 
fa ^~4»l<Ti (C.14) 



C.4. Statistical error on the estimated differential pistons 

To compare the various schemes, we are interested by the error 
on the differential piston x, ; = jCj - Xj, which corresponds to the 
error on the correction applied to the delay lines: 

B 

x i = J]\V ik fa (C15) 

k=\ 

Given the definition of fa (Eq. IC.14l ), the statistical error on these 
terms is cr(fa) = 1. We finally get the quadratic error <x 2 . on the 
corrected differential piston: 
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(C.16) 



